Astronomy & Astrophysics manuscript no. paper 


© ESO 2009 


June 4, 2009 





A dam around the Water Fountain Nebula? 

The dust shell of I RAS1 6342-381 4 spatially resolved with VISIR/VLT 

T. Verhoelsti'2 *, L.B.F.M. Waters'-'^, A. VerhoefF*, C. Dijkstra^ H. van Winckeli , J.W. Pel^ and R.F. Peletier^ 

' Insituut voor Sterrenkunde, K.U. Leuven, Celestijnenlaan 200D, B-3001 Leuven, Belgium 

e-mail: ti j 1 . verhoelst@ster . kuleuven . be 
' University of Manchester, Jodrell Bank Centre for Astrophysics, Manchester, M13 9PL, U.K. 

Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA 
* Astronomical Institute 'Anton Pannekoek", University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands 
^ Kapteyn Astronomical Institute, Landleven 12, 9747 AD Groningen, The Netherlands 



o 
o 

(N 



0^ 

6 



o 

O 

o 

ON 

1. Introduction mate the fast Post-AG B wind have now been observed around bi- 

J> nary Post-AGB stars dPeroo et al. 2006, 2007) and observations 

. ^ Two different mechanisms for the shaping of bipolar and multi- ^^e dust at the heart of the Ant Nebula ( Chesneau et alfcOOTh 

polar (Proto-)Planetary Nebulae (PNe) are cuiTendy under in- consistent with a disk interpretation. On the other hand, op- 

H tense investigation (e.g., |Meixner et alj | 2004i) . The first is in ^j^^j ^^^.j^ observatio ns of some P PNe show evidence for 

. , con£e£t based on the Generalized Interacting Winds model sharping bv (precessing) iets (ISahai et al.ll2007l) . The presence of 

(lBalic^Ll987) and suggests the colhmation of the fast, spheri- ^ ^^^^^ ^^^^^ ^j^^ suggested in these objects, but it is hitherto 

cal wind of the hot central star by an older, slow, axi-symmetnc ^^^^j^^ ^^^^^^^ ^^^^^ ^^^i ^^^^^^^ ^^^^^ observed with 

and possibly warped, AGE (super-)wind (e.g. Icke 2003). The jy^jj^ interferometry 
axi-symmetric shaping of the earlier AGB wind may require the 
central star to be in a binary system. 

The other scenario postulates that the lobes are in fact cav- IRAS 16342-38 14, hereafter IRAS 16342 is a young PPN. In 

ities blown in a preceding spherical AGB wind by (precessing) spite of its exti-emely red SED, wit h crystalline silicate features 

jets dSahai et al.ll2005. 2007,) . The jets are then most likely the in absorption up to almost 45 /urn jPiikstra et al.ll2003h . optical 

consequence of (re-)accretion of circumstellar material, maybe HST images s how a bipolar re flection nebula with a dark equa- 

also in a binary system. This scenario is especially attractive to torial waist (Sahai et al1ll999 ', hereafter STM99). According to 

explain very young PPNe, in which the central star is not yet hot STM99, the lobes are cavities blown by a bipolar molecular 

enough to exhibit a strong line-driven wind. Both hypotheses are jet, observed as high-velocity lines of water (hence the designa - 

not mutually exclusive. tion "water fountain nebula") and OH dLikkel & Morrislll988h . 

With the advent of powerful mid-IR interferometers, espe- within the remains of a low mass- loss AGB wind w hich pre- 

cially MIDI at the VLTI, the type of dust disks which could colli- ceeded the recent superwind phase. Sahai et aTl ( 12005 ^) interpret 

local density enhancements in the bipolar cavities as due to pre- 

Send offprint requests to: T. Verhoelst cession of the jet. The dark waist, origin of the bulk of the en- 
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ABSTRACT 

Context. Bipolar morphologies in Planetary Nebulae (PNe) are believed to be closely linked to binary central stars. Either by coUi- 
mating a fast stellar wind or by driving a jet via accretion in the central system, dusty torii or stable disks may be crucial ingredients 
for the shaping of PNe. 

Aims. We study the dust distribution in the very young Proto-Planetary Nebule (PPN) IRAS16342-3814, also known as the Water 
Fountain Nebula, which is known to show strong bipolar characteristics in the shape of two reflection lobes, and high-velocity coUi- 
mated molecular outlfows. 

Mefhods. We use the new Mid-IR (MIR) instrument VISIR on the Very Large Telescope (VLT) both in imaging and spectroscopy 
mode at wavelengths from 8 to 13 pm. 

Results. We present the first spatially resolved MIR observations of a dusty evolved star obtained with VISIR and find that the im- 
proved spatial resolution contradicts previous claims of an elliptical brightness distribution at the heart of IRAS 16342: we find the 
waist region to be dark even in the MIR. We show that the filling angle of the obscuring dust lane, which is made mostly of amorphous 
silicates, is very large, possibly even close to a spherically symmetric superwind as seen in OH/IR stars. 

Conclusions. We conclude that, in contrast to the multitude of recent dusty-disk detections in Post-AGB stars and PNe, IRAS 16342 
does not show this extreme equatorial density enhancement, at least not on the scale of the dusty environment which lends the object 
its IR appearance. Rather, it appears that the observed precessing jets are shaping the bipolar nature in the remains of a spherically 
symmetric AGB superwind. 

Key words. Techniques: spectroscopic - stars: AGB and post- AGB - stars: circumstellar matter - stars: individuaklRAS 16342-38 14 
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which completely obscures the central star in the line-of-sight to 
the observer. 

iDiikstra et all (|2003|) . hereafter DKW03, present infrai-ed 
ISAAC and TIMMI2 images from 3.8 to 20 yum and find only 
at the shortest wavelengths a bipolar nature such as in the HST 
images. At longer wav elengths, the s ource appears elliptical, as 
was already found by Meix ner et al.l ([199^. 

In this paper, we present a new MIR image and an N-band 
spectrum of IRAS 16342 obtained with the VISIR instrument on 
the VLT, both of which resolve the dust structure and force us to 
question our understanding of this (type of) source. 



2. Observations 

The observations presented here were obtained as part of the 
Dutch GTO on VISIR, the MIR imager and spectrometer in- 
stalled at the Cassegrain focus of Melipal at the VLT. An image 
was made in the SiC filter (Ac = 11.85yum, A/l - 2.34jum) with 
a 0.127 ai-csec PFOV in the night of 20 March 2006, bracketed 
by 2 PSF observations (HD 146051, M0.5III). Low resolution 
spectroscopic observations in 4 bands (centered at 8.5, 9.8, 1 1.4 
and 12.2 jum) were made in the night of 17 March 2006 with a 
32x0.75 arcsec slit oriented 55.4° East of North, and a 0.127 arc- 
sec PFOV. The calibrator was again HD 146051. science obser- 
vations ranges from 1.03 to 1.05 and that standards from 1.08 to 
1.12. The airmass of the imaging ranges from 1.15 to 1.25. The 
visual seeing was about 0.60 arcsec for the imaging observations 
and 0.90 arcsec at the time of the spectroscopy. 

For the data reduction, we use the pipeline offered by ESO, 
version 1.3.7. The synthetic N band spectrum of the calibrator 
is co mputed using a marcs atmosphere model (Gus tafsson et alj 
11975 , and further updates ), including SiO op acity. The stellar 
parameters are taken from iBorde et al.l ( |2002|) . The telluric cor- 
rection of the N band spectra is performed by solving for in- 
strumental efficiency and atmospheric optical depth using the 2 
observations of HD 146051. The resulting numbers are then in- 
terpolated in airmass to match the conditions at the time of the 
science observation. As absolute calibration, we scale the VISIR 
spectrum in such a way that the simulated inband power over the 
IRAS 12 jum filter matches the observed value of 16.2 Jy (see 
Sect. |3.1| for a discussion on possible slit losses). 



3. Bipolar even in the IVIid-IR 

The image deconvolutio n is p e rformed wit h both a maximum- 
likeUhood (Richardso nI II972I: iLucvl 1 1 974 and a maximum- 
entropy method ( Agm on et alj|1979l) . Differences between the 
results of both methods are minimal. The raw and deconvolved 
images are shown in Fig.[T] 

We observe a double-peaked intensity distribution with a 
separation of 0.92 arcsec and a Position Angle (PA) of 66° East 
of North. Component W (West) contains about 3/4 of the total 
flux, component E (East) the remaining quarter TIMMI2 im- 
ages in the N and Q bands observed by DKW03 do not resolve 
2 separate emission peaks, but show instead an elliptical shape 
with a major axis oriented more or less like the separation vec- 
tor observed at shorter wavelengths. The PA of the object in our 
VISIR image is compatible with the values found in the optical 
and near-IR, but the separation we find at 11.85 fj.m is slightly 
smaller than that at shorter wavelengths. We conclude that the 
elliptical intensity distribution observed by Meixner et al. ( 1999) 
and DWK03 is due to insufficient spatial resolution in their ob- 
servations. 





Fig. 1. Upper panel: The raw image of the Water Fountain 
Nebula in the SiC filter (11.85 yum). The colour scale is linear 
The direction of the bipolar optical reflection nebula is also in- 
dicated (solid line, from STM99). The dotted lines represent the 
position and width of the slit used for the spectroscopic observa- 
tions. Lower panel: The same image but deconvolved using the 
PSF observation of HD 14605 1 and with the colour scale follow- 
ing I*'^^\ The solid and dotted lines have the same meaning as 
in the upper panel. 



In the deconvolved image, it is apparent that both lobes do 
not point at a common centre of symmetry. This was also ob- 
served by Sahai et al. ( 1999) in the optical HST images. 

Using the central starparameter£|of DKW03 (T^s = 2670 K 
and Rt - 372 Rq), the distance of 2 kpc and inclination / = 40° 
from STM99, and assuming that the dust temperature scales 
with the distance r from the central star as = reff(2r/7?,)"°'* 
(e.g. .Herman et al. 1986), we find that ~ 160 K at 900 AU, 
i.e. the distance between the central star and the approximate 
center of either lobe. Both numbers are uncertain to at least 20%, 
but what follows does not depend so sensitively on these num- 
bers. Using Wien's law, we find a peak of the energy distribution 



' These parameters correspond to a star at the tip of the AGB, but 
we can not rule out a more evolved central star. This would imply a 
slightly higher temperature for the small dust grains, but has no major 
consequences for the results discussed here. 
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Fig. 2. The reduced and calibrated, but not yet collapsed, spectra 
(4 bands) pasted together and normalized to the continuum level 
as derived from a spline fit to the ISO-SWS spectrum (DKW03). 
The colour scale is linear Perusal of this figure already reveals 
that the dark waist reduces the flux level of both lobes around 
lOfim, but more so for the Eastern lobe. The solid white lines 
represent contours of constant opacity following the model pre- 
sented in Sect.m 



just to the red of the N-band. Given the spectral shape of the 
lobes over the N-band presented in Sect. |3.U we conclude that 
the N-band image and spectrum are dominated by thermal emis- 
sion from the dust in the bipolar lobes, which is heated directly 
by the central star. The material in between, which is obscuring 
the central star, remains dark even in the MIR, indicating that 
along the line-of-sight, optical depth r = 1 is reached already in 
the cold (r < 100 K) outer regions of the dusty environment. 

3. 1 . The spectra 

The four reduced 2D spectra (at 8.5, 9.8, 11.4 and 12.2 fim) are 
shown in Fig.|2] one above the other following increasing wave- 
length. As shown in Fig.[Tl slit losses should at most be of the 
order of a few percent. From this 2D spectrum, we can derive the 
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Fig. 3. In black, we show the field-integrated pipeline-extracted 
spectrum calibrated as discussed in Sect.|2l For comparison, the 
ISO-SWS spectrum and spline-fit continuum of DKW2003 are 
shown in grey and dotted lines respectively. Shown in red and 
green are the separated spectra of the Western and Eastern lobes 
respectively. The grey box indicates the wavelength region of 
strong atmospheric ozone absorption. The calibration procedure 
presented in Sect.|2]resulted in reliable data even at these wave- 
lengths. 



total field-integrated spectrum, shown as the black solid line in 
Fig. [3] which is of higher quality than the ISO-SWS spectrum of 
DKW03 (shown in grey, scaled to match the absolute IRAS flux 
level) but has the same global shape and no particular spectral 
feature^ in contrast to what is observed at longer wavelengths, 
we don't find a clear indication of crystalline silicates (e.g. the 
Forsterite feature at 1 1.3 fim) within the N band 

More interesting is the possibility to separate the Eastern and 
Western lobes in this 2D spectrum. The slit PA of 55.4° is very 
close to the object PA of 66° derived from the VISIR image pre- 
sented in Sect. [3] To extract position, width and flux levels for 
the individual spectra, we performed a row-wise fit with a dou- 
ble Gaussian intensity distribution. The 2 lobes can be fully sep- 
arated through the entire wavelength range. 

The resulting wavelength-dependent positions and widths 
for both lobes are shown in Fig.|4] The wavelength-dependence 
of the separation between both lobes clearly resembles the opac- 
ity profile of amorphous silicates, the primary consituent of the 
dust grains surrounding the central star. The presence of such a 
wavelength dependence shows that the extinction is not due to a 
uniform screen of dust (ISM or detached relic of an early-AGB 
wind) but instead that, for any given wavelength, the opacity de- 
creases with distance from some point in between both lobes. 
This is also confirmed by the detection of a varying width in the 
marginally resolved individual lobefl Since both lobes are af- 
fected by the dark waist, we find a disagreement with the model 
presented by DKW03, in which the flat dusty disk does not result 
in any obscuration of the Western lobe. 

The individual spectra of the Eastern and Western lobes are 
presented in Fig. [3] and in Fig.|5]we display the ratio between 
both. The individual spectra show that both lobes emit an equally 
cold spectrum, and both exhibit the silicate feature in absorption. 
This is again in disagreement with the low filling angle of the 
model by DKW03 which predicts a dust emission spectrum for 
the Western lobe. 

The flux ratio shows that the Eastern lobe suff'ers stronger 
attenuation than its Western counterpart. Again, the major fea- 



- The absorption line at 12.5 - 12.6/im is telluric. 

' The standard star observations show a width of 0.27 : 



: 0.03 arcsec. 
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Fig. 4. The dotted and dashed Hnes represent the position of the 
E and W lobes respectively, minus the position of their contin- 
uum at 8.5 jum and offset by 1.5 arcsec for clarity. The solid line 
represents the separation between component E and W as a func- 
tion of wavelength. The black diamond at 1 1 .85 fim indicates the 
separation measured in the 2D image presented in Sect. [3] The 
two curves at the bottom show the width of the individual com- 
ponents . 
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Fig. 6. The toy model used to estimate the filling angle of the 
obscuring dust structure. The density in the dust structure drops 
with r^, where r is the distance to the central star. The inclination 
is taken from SAH99. (IS,y) is the intersection between a line- 
of-sight and the inner edge of the dust structure. The lobes emit 
thermal radiation at 160K. 
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Fig. 5. The observed flux ratio between component E and W 
(black solid line) and the model presented in Sect.|4](red solid 
line). Also shown are typical absorption (dotted lines) and scat- 
tering (dashed lines) cross sections for iron grains, and for a mix 
of small and large silicate grains. The scattering by iron grains is 
negligible at these wavelengths. The red dash-dot line represents 
tiie /c(/l) profile of the best-fit dust composition. 



ture is that of amorphous silicates. A proper identification of the 
type of silicates would also require the observation of the 18 [im 
band. The flux ratio Fy^^lFy;sN ~ 0.4 (instead of unity) at the 
"continuum" wavelengths of 8.5 and 13 fim indicates achromatic 
attenuation, which is most likely due to metallic Fe, a crucial 
ingredient in the SED modellin g of O H/IR stars and obscured 
RSG jKemper et alJl2002tlHarwit et al.ir200lh . The observed ab- 
sorption towards 8 //m can be attributed to amorphous silicates, 
but only through scattering by fairly large grains. This suggests 
that a fraction of the obscuring dust is contained in micron- 
sized grains, which was already suggested for the dark waist by 
DKW03, albeit based on a geometrically inadequate model as 
discussed above. 



Having found that the obscuring dust structure must have a 
significant filUng angle, we try to quantify this in the following 
section. 



4. Interpretation: filling angle 

Located at a similar distance as many known Post-AGB bina- 
ries, i.e. a few kpc, IRAS 16342 appears much larger in the mid- 
IR than the stable dusty disks around those stars: 1 arcsec vs. 
40 mas respectively. Equivalently, the observed strong extinc- 
tion toward the Eastern lobe requires a significant dust column 
at least 500 AU above the midplane, and this is not possible in a 
stable Keplerian dust disk with a reasonable mass. The dust must 
therefore have some radial velocity. 

With a simple geometrical model (see Fig.|6]l, which is not 
intended to be a detailed representation of the source with all its 
peculiar characteristics, it is possible to get an estimate of the 
filling angle and dust mass of the obscuring dust structure from 
the spatially resolved N-band spectrum. The mid-IR data pre- 
sented here have as advantage over the optical and near-IR data 
presented hitherto that the modelling is much less sensitive to 
scattering. Including scattering in the analysis would add a new 
level of complexity (e.g. the central source spectrum, the scat- 
tering properties of the dust at short wavelengths, the scattering 
phase function), and is not required to constrain the basic struc- 
tural parameters. We leave the construction of a full "colour" 
image to a future paper 

In the assumption of an outflow with a biconical cavity, the 
optical depth r along a line-of-sight at a distance y from the cen- 
ter is calculated as: 



T(y,A) 



E,W 



f 

Jis 



dx 



)is ^ 
Pqk(A) 

y2 
y2 



-r 



[arctan(x/y)] 



>'tan(f-f±i) 
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It depends on (1) the filling angle 6 of the dust structure, (2) 
the dust density po at unit radiufl and (3) the dust mass frac- 
tion caught in small and large silicate grain £1 (0.1 and 5yum), 
and in metallic iron, implicit in the wavelength-dependent cross 
section k. The "background" radiation by the lobes is assumed 
to be thermal emission by small dust grains at 160K, as derived 
in Sect.[3] Whether these small dust grains fill the cavities or 
form the boundary with the denser dust structure cannot be de- 
rived from our current observations. We do not assume an actual 
density or temperature distribution for the individual lobes: the 
lobes appear marginally resolved in our observations, but insuffi- 
ciently so to derive such a distribution. The comparison between 
model and observations is therefor done at the "nominal" dis- 
tances of yo = 0.45 arcsec (900 AU) from the center of the ob- 
ject. We determine the unknowns by requiring reproduction of 

(1) the very red SED up to 20 /im with total obscuration of the 
central star and hot inner dust (from the optical to the mid-IR), 

(2) the wavelength-dependent flux ratio between both lobes, and 

(3) the moderate observed extinction towards the Western lobe 
(tio/jiii ~ 0.5, derived from the depth of the silicate absorption 
feature in the Western spectrum). The wavelength dependence 
of the dust opacity is mostly constrained by (2). The dust den- 
sity is determined by (1) and (3). 

We find that a filling angle of at least llOdeg is required 
to produce observable extinction towards the Western lobe. The 
best fit is obtained with a 145 deg filling angle. With a dust 
species mass ratio of - 5.5, a silicate size dust mass ratio 

of = 1.5, and a mass-loss rate of M ~ 1 x lO^^Moyr"', 

we find good agreement with the N-band observations (Fig.|5]l 
and lack of flux at shorter wavelengths (toward the central star, 
our model predicts r > 100 up to 30/im). However, we can not 
exclude an even larger filling angle. 
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5. Discussion and conclusions 

We find that, in spite of previous classifications as an elliptical 
MIR source, IRAS 16342 in fact appears as a bipolar source up 
to at least 13/^m, if studied with sufficient angular resolution. 
Almost all of the N-band flux originates in these 2 lobes, with 
an apparant temperature of roughly 160K. The dusty "waist" 
region remains dark even at these wavelengths. From the vary- 
ing extinction in front of both lobes as a function of wave- 
length, we derive that the dust structure has a large filling angle 
(~145°), much larger than the stable Keplerian dust disks re- 
cently seen around other binary Post-AG B stars (typically 40°, 
Ide Ruvter et~aLl 120061; iDeroo et al.l l2007l) , and this in spite of 
the strong indications that also IRAS 1 6342 has a binary cen- 
tral system. Instead, it resembles more the spherically symmet- 
ric superwind of an OH/IR star, but with cavities blown by jets. 
Unfortunately, the dam around this fountain is too high to see 
the engine driving the bipolar jets. It will be interesting to see 
whether all members of the class of Water Fountain Nebulae dis- 
play this kind of dust structure. 



Acknowledgements. The authors would like to thank the anonymous referee 
for many valuable comments, and E. Lagadec for a careful reading of the 
manuscript. 



The density follows a r ^ distribution, with a 15 km s ' outflow ve- 
locity. 

^ We use only these 2 discrete sizes since there is no reason to believe 
that a classical power law size distribution is valid for this dust structure. 



